Particle escapes in an open quantum network via multiple leads 
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Quantum escapes of a particle from an end of a one-dimensional finite region to N number of semi-infinite 
leads are discussed by a scattering theoretical approach. Depending on a potential barrier amplitude at the 
junction, the probability P{t) for a particle to remain in the finite region at time t shows two deferent decay 
behaviors after a long time; one is proportional to N'^ /t^ and another is proportional to l/{N'^t). In addition, 
the velocity V{t) for a particle to leave from the finite region, defined from a probability current of the particle 
position, decays in power ~ 1/t asymptotically in time, independently of the number A'' of leads and the initial 
wave function, etc. For a finite time, the probability P{t) decays exponentially in time with a smaller decay rate 
for more number A'^ of leads, and the velocity V{t) shows a time-oscillation whose amplitude is larger for more 
number of leads. Particle escapes from the both ends of a finite region to multiple leads are also discussed by 
using a different boundary condition. 

PACS numbers: 05.30.-d, 03.65.Nk, 73.23.-b 



I. INTRODUCTION 

The particle escape is a typical nonequilibrium phe- 
nomenon in open systems. It is a cuirent of particles from 
a region where the particles were initially confined. The con- 
cept of escape has been used to describe a variety of physical 
phenomena, such as a-decaying nucleus iQl-lst], chemical reac- 
tions as Kramers' escape problem 10-05 etc. Some dynamical 
properties, like chaos ||7H9[], the recurrence tim e jlOtl . the first- 
passage time Sl^, transport coefficients ifTlUlSll . etc.. have 
also been investigated via escape behaviors of particles. The 
particle escapes were discussed in many types of systems, for 
example, stochastic systems J|4j6jl, classical billiard systems 
|7,[M 14 jl^ . map systems lUlSESi and quantum systems 



Escape phenomena in open systems cause a decay of var- 
ious quantities by a particle current from a finite region. For 
instance, the probability for particles to remain in the initially 
confined region, which we call the survival probability in this 
paper, would decay in time, if particles continue to be flowed 
out from the initial region. Such decay properties in escape 
systems have led to some interesting results and conjectures. 
For example, in classical mechanical systems with a particle 
escape via a small hole, the survival probability would de- 
cay exponentially in time if the dynamics is chaotic, while 
it would decay in power for non-chaotic systems 
On the other hand, in many quantum mechanical systems, 
the survival probability decays in a power asymptotically in 
tinie 0l[T8ll20l - l23ll with an exponential decay for a finite time 
lUmil^K and values of the power in the d ecay vary by initial 
conditions ll24ll or particle interactions Ii25.i26ll . etc. 

In this paper we discuss particle escapes in quantum me- 
chanical networks as an example of open dynamical sys- 
tems. The quantum network system is also called the "quan- 
tum graph," and is constructed by connecting finite and infi- 
nite narrow wires like a network, and have been widely used 
as models to describe mesoscopic transports like Aharonov- 
Bohm type of effects ll27ll28ll . resonance tunnellings ll29l [soll . 
cuiTent splitters 131143311 . chaos and diffusion 134 13511 . etc. 



Steady electric currents in open quantum network systems are 
described by quantum scattering theory ll36l - [39ll . This kind of 
quantum systems with narrow wires could be experimentally 
realized as a combination of atomic or molecular wires or as 
a graph-like structure on the surface of a semiconductor by a 
recent development of nano-technology ll40t - l42ll . 

An important feature of network systems is an effect of 
current splitter at a network junction. In order to consider 
such a splitting effect of currents in quantum escapes as sim- 
ple and concrete as possible, we consider particle escapes 
from a finite one-dimensional region with the length L via 
A^-number of semi-infinite one-dimensional leads. The mul- 
tiple leads is connected at one end of the finite region with 
a potential barrier amplitude A at the junction, and we im- 
pose the fixed boundary condition at another terminated end 
of the finite region. As a theoretical approach to describe 
particle escapes in such a quantum network, we use a quan- 
tum scattering theoretical approach, by which we consider 
concretely decay properties of two quantities to characterize 
the quantum escapes. The first quantity of decaying in the 
particle escape is the survival probability P{t) for the par- 
ticle to remain in the finite region at time t. We show an- 
alytically that the survival probability P{t) depends on the 
number A^ of attached semi-infinite leads as the A^'^ pro- 
portional ratio, i.e. Wmt^j^oo P{i) /[P{t)\N=i\ = N"^, in 
the case of A 7^ — fi,^/(2mL) (m: the mass of particle, 
2Tih: the Planck constant). This means that after a long 
time the particle can remain more inside the original finite re- 
gion by connecting more semi-infinite leads. Besides, in this 
case the survival probability P{t) decays in power ~ \/t^ 
asymptotically in time. In contrast, in the case of A = 
—h'^/{2mL) we obtain the A^^ inversely proportional ratio 
lm-it^+rx,Pit)/[P{t)\N=i] = 1/A^2 foj. the survival proba- 
bility P{t), meaning that after a long time a particle escapes 
more by connecting more leads. Here, the survival probabil- 
ity P{t) decays in power ^ 1/t asymptotically in time, dif- 
ferently from the case of A 7^ —h'^/{2mL). We also discuss 
finite time properties of the survival probability numerically, 
and show that for a finite time the survival probability decays 
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exponentially for a longer time with a smaller decay rate by 
connecting more number N of leads. As the second quantity 
to investigate decay properties in quantum particle escapes, 
we consider the velocity V{t) for the particle to escape from 
the finite region, which is introduced by the equation of con- 
tinuity for the particle position probability. We show analyt- 
ically that the escape velocity V{t) behaves asymptotically 

in time as V{t) L/t for A 7^ -ftV(2mL) and 

V{t) L/{M) for A = -ff/{2mL), which are inde- 

pendent of the number N of leads and the initial wave func- 
tion, etc. It is also shown by numerical calculations that for 
a finite time the escape velocity V{t) oscillates in time, and 
takes a larger magnitude of the oscillations for more num- 
ber N of attached leads. Furthermore, by using a different 
boundary condition, our scattering approach to quantum net- 
work systems also allow us to discuss particle escapes from a 
finite region with the length 2L whose both end are connected 
to N number of semi-infinite one-dimensional leads. 

The outline of this paper is as follows. In Sect. [Ill the quan- 
tum network system with a finite wire connected to multiple 
leads is introduced. Based on the equation of continuity for 
particle position probability, we introduce a probability cur- 
rent whose conservation imposes boundary conditions at the 
junction of leads and at the terminated end of the finite wire. 
These boundary conditions specify the scattering states of this 
system, from which we describe the time-evolution of a wave 
function of the system. In Sect, [nil we introduce the sur- 
vival probability P{t) from the particle position probability 
and the escape velocity V{t) from the probability current in 
the quantum network system, and discuss these decay proper- 
ties. Finally, we give some conclusion and remarks in Sect. 

m 



II. QUANTUM SCATTERING APPROACH TO NETWORK 
SYSTEMS WITH MULTIPLE LEADS 

A. Quantum network system and boundary conditions 

We consider quantum network systems consisting of a finite 
one-dimensional segment with a length L whose end is con- 
nected to N number of semi-infinite one-dimensional leads. 
(See Fig. [T]as a schematic illustration of this network.) We 
call the finite segment with the the length L the region 
and also call the j-th semi-infinite segment of lead the region 
(n 1, 2, • • • , N). In each region we put the cc'"'- 
axis of coordinates with the origin O at the junction of leads, 
in which the positive direction of the a;'") -axis is taken from 
the origin O to the region (n = 0, 1, 2, • • • , N). 

We introduce the wave function ^'("^(a:, t) of a particle in 
this quantum network at time t and the position x G /C'"' 
for X > 0. From Schrodinger equation for the wave func- 
tion 4''^"' (x, t) with a real potential we derive the equation of 
continuity for the particle position probability density 
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FIG. 1: A quantum network with a junction for multiple leads. It 
consists of a finite region £''^' with a length L and the semi-infinite 
one-dimensional region (n = 1,2, ■ ■ ■ , N) connected to an end 
of the region C^'^^ . The a::^"'-axis of coordinates is taken in the region 
with the origin O at the junction and its coordinate satisfies 
x'"' > (n = 0, 1, 2, • ■ • , iV). The coefficients a|"'^j (^^(",1)) 
the amplitude of the incident wave to (the scattered wave from) the 
junction in the region by a plain wave injected from the j-th 
lead with the wave number k. 



dt dx 

in which the local velocity i;*^"^ [x, t) at the position x G 
and the time t is introduced as 



{x, t) = — Im 



dx^ 



(3) 



with the mass m of particle and the Planck constant 2-Kh ||43|] . 
Here, Im[Ar] means the imaginary part of X for any complex 
number X. 

For describing the quantum state at the junction O, we im- 
pose that the wave function of the system is continuous at any 
position, including at the origin O so that the boundary condi- 
tions 



^'(0,t) = lim ¥°\x,t) = lim ^'(^'(a;,^ 
>+0 3;— >+o 

= ••• = lim ^^^\x,t) 
2;— y+O 



(4) 



are satisfied at any time t. We further assume that there 
is no net particle current source at the junction O, namely 

lim,^+o T.n=o P^"' ^)«^"^ = 0- leading to 



(1) 



N 

lim Vu(")(a;,i) = 

n=0 



(5) 
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for the local velocity (|3]l, noting that by the condition (|4|l the 
position probability density lim 
dent of the region number n 
Eq. ^ is rewritten as 



+0 p''"'' {x, t) is indepen- 
Using Eqs. Q and (01, 



lim 



N 

E 



d¥'^\x,t) 
dx 



(6) 



with a real constant A Bo, l39l 14511 . It is important to note 
that for the case of iV = 1 the condition (|6]l is equivalent to 
the boundary condition for a delta functional potential with 
the amplitude ft^A/(2m) in an one-dimensional one-particle 
system. In this sense, we regard the real constant 



A = 



2m 



(7) 



with the constant A appearing in Eq. (|6]l as a potential barrier 
amplitude at the junction O. Finally, we impose that there is 
no particle current source at the end .t'^^^ = L of the region 
£(0), so 



,,(0) 



(L,t) = 0, 



namely 



dx 



(8) 



(9) 



with a real constant /^j, similarly to Eq. ^ 144 14511 . 



B. Scattering states and matrix 



suppressed the suffix j for the scattering matrix Sk, because 
as shown later in Eq. ( fT2] i the scattering matrix is indepen- 
dent of the suffix j. The energy eigenstate (fTOl i. which has 



nonzero value even in the infinite region x 



in) 



.of/:("). 



n = 1, 2, • • • , A^, is the scattering state of the quantum net- 
work system. 

The scattering state ( fTOl i as a special example of the wave 
function must satisfy the conditions ^ and 

leading to the specific form of the scattering matrix Sk as 



N + l + ii 



-U -I. 



(12) 



Here, / is the (A^ + 1) x (A^ + 1) identical matrix, and U is 
the {N + 1) X {N + 1) matrix whose all elements are given 
by 1. [See Appendix lAl for an derivation of Eq. ( fTTT i with the 
scattering matrix ( fT2l i. ] Noting the relation = {N + 1)U, 
the scattering matrix (fTSl i is shown to be an unitary matrix 
satisfying the relation 



s}^Sj 



SkSl = / 



(13) 



with the Hermitian matrix S], of Sk, supporting the relation 



a{j,k)\^ by Eq. O- 
On the other hand, the condition ^ leads to the relation 



1 



21*1 ^2i/cL. (0) 



(14) 



We assume that the quantum network introduced in the pre- 
vious subsection III Al consists of ideal leads, i.e. there is no 
potential in any of the region ?7, 0,1,--- ,iV except 
at the junction. In this case, the energy eigenstate $|j^^-j(a;) 

at the point x in is presented by a superposition of plain 
waves with the wave number k (> 0) as 



between the amplitudes a^^\-f and ^ for the scattering state 
in the finite region C^^\ The condition ( fT4l i guarantees the 



same magnitude 



^ U k) 1^ plain waves prop- 

agating to the opposite directions with each other in the finite 
region The relation (fT4l i includes two special cases: 



(n) 



(x) 



(ri - 

a; . , , e 



ikx 



6^") e 



ikx 



(10) 



n = 0,1, - ■ ■ ,N, corresponding to the energy eigenvalue 
Ek = h^k^ /{2m). Here, we introduced the suffix j in 

ferent states with the same energy Ek as discussed later in 



the quantities <i>^J^j,^ (x), a^^-'k) and to distinguish dif- 



detail. The constants ci-^fk) ^{j i-) '^l' 
garded as the amplitude of incident wave to and scattered 
wave from the junction O, respectively. The {N + 1) dimen- 

is 



sional column vector bi 



■■■b 



'U,k) - y^U,krijM) "U^k) 

connected to the {N + 1) dimensional column vector a(j k) 



^ U,k)"-U,k) 



•«0-,fc) ) as 



b 



U.k) 



Ska 



kO-U,k) 



(11) 



with the scattering matrix 5*^. Here, the notation de- 
notes the transpose of X for any vector X. In Eq. ( fTTT i we 



(0) 

{3,k) 

(0) 



2ikLA0) 
''U-k) 

2ikL,{0) 
°U,k) 



for fi —>■ +00, (15) 
for/i = 0, (16) 



meaning that the limit ^ — > +oo imposes the fixed (Dirich- 
let) boundary condition (l>|°''j,^ (L) = for scattering states 
and the case of /i = imposes the open (Neuman) boundary 
condition at the end x'^'^^ = L of the finite region C^^\ It 
is important to note that the open boundary case /i = and 
the fixed boundary case /i +oo can also describe states in 
the quantum network consisting of a finite lead with a length 
2L whose both sides are connected to N number of multi- 
ple semi-infinite leads, as shown in Fig. [J] Here, the open 
boundary case i^i = (the fixed boundary case — >■ +oo) cor- 
responds to the quantum state described by the wave function 
with the inversion symmetry (anti-symmetry) at the center of 
the finite region with the length 2L in such a system. 
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FIG. 2: A quantum network with two junctions for multiple leads at 
the both ends of the finite region with a length 2L. 



By Eqs. (fTTT i. (fT2l i and ( fT4l i the coefficient &["^ is given by 



.,(0) 



„("') 



iV + 1 + i| + (A^ - 1 + ^^e^^'^'f' ' 



(17) 



2(1 + i±i|e2''^-^ 



Z^n' = l "(j.fe) 



(18) 



for ?7 = 1, • • • , iV. Eqs. (O, ([TtIi and ([TSll shows that the co- 
efficients al?,,, and 6'"'^.,, n = 0, • • • , are given from the 



(n) 

coefficients a^^ n = 1, • • • , A^. We can still take any value 



for a 



1 , • • • , A^, as the A^ number of incident wave 



amplitudes from infinite regions a:'"^ — >■ +oo, n= 1, • • • , A^, 
and this voluntariness leads to necessity of the suffix j to dis- 
tinguish degeneracy of the energy eigenstates ^(J^-) (2;) corre- 



sponding to the energy E^. In this paper, we introduce a^^ 

as 



(19) 



for n = 1, 2, • • • , A^ with a constant ^j. Eq. ( fT9] l means that 
the state (a^) is the scattering state produced by the inci- 
dent wave from the infinite region x^^^ — +00 of only. 
We define the inner products 



[X\Y) = lim 

e->+0 



dxx(°)(.T)*y(°)(x) 
+ lim / A:(")(.T)*y(")(x) 

'^^'^^ n=l •''^ 
N 

= y / dx x'^"\xyY^''\x) (20) 

for any quantum states X and Y, whose values at the posi- 
tion X in the region are given by X^'^\x) and y^")(a;), 
respectively. Here, X* with the asterisk denotes the complex 



conjugate of X for any complex number X. For this inner 
product (I20I 1 it is essential to note that the scattering state 
$|"^^.j(a;) satisfies the orthogonal relation 

i'^U.k) I 'Pu'M)) = ^n'S{k-k') (21) 
for fcfc' > 0. Here, we specified the constant in Eq. ( fT9l ) as 

1 



(22) 



independently with respect to the suffix j = 1, 2, • • • ,N, so 
that the coefficient of the term Sjj'6{k — k') in the right-hand 
side of Eq. ( 1211 1 becomes 1. The proof of Eq. ( I2TI ) is given in 
AppendixlB] 



C. Green function and the time-evolution of wave function 

Now, we assume that the initial quantum state ^I* (0) at t = 
is expanded by the scattering states, as 

N „ + oo 

¥-\x,0)^Y. (23) 

at the position x in with a constant By Eqs. ( l20b . 

( I2T] ) and ( I23] ) the constant ^(j.fe) is expressed as 



%^) = (%^) I *(o)> 



(24) 



for k > 0. The initial condition ( |23] l implies that the initial 
state ^'(0) does not include bound states with discretized en- 
ergy spectrum which are orthogonal to the scattering states 

The wave function of the system at time t > is given by 
applying the time-evolutional operator cxp{—iHt/h) with the 
Hamiltonian operator H to the initial state ( |23] l. and we obtain 



(a;, = y / dkA 



-iEkt/h^{n) 



ix) 



3 = 1 

N 



= y / dx' G("'"')(x,a;';t)4'("'H:E',0) 

(25) 

by Eqs. m. (El and i?$S"i)(a;) = Ek^^[^)^){x). Here, 
; <) is the Green function defined as 

G("'"')(a;,a;';t) 



N 

E 



dk $ 



S)(-)<j';;y)(-')*e 



("') (rr'Yp-^Ekt/h 



(26) 



In other words, we can calculate the Green function ( |26] ) via 
the scattering state ([TOl) with the coefficients (O, ([T7]i, (fTsT i 
and (fT9] l. leading to the wave function at time t by Eq. dZSl l. 
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III. ESCAPE BEHAVIOR OF A PARTICLE VIA 
MULTIPLE LEADS 

A. Quantum escapes in a network system 

In this section we consider behaviors of a particle to escape 
from the finite region to the semi-infinite regions C^^', 
n = 1, 2, • • • , A^. In order to discuss such phenomena, we set 
the initial wave function at t = to take nonzero value only 
in the finite region 

'J'(")(.T,0) = forn = I,--- ,iV. (27) 

By Eqs. (l25l l and (IZTl i the wave function ^I/C) (x, i) in the 
region at time t is expressed as 

¥°\x,t)^ [ dx' G^°^°\x,x';t)¥°\x\0) (28) 
Jo 

I 



by using the functions G'^°'°^x,x'; t) and vI/(o)(aT', 0) for this 
finite region only. 

As shown in Appendix |C] the Green function 

G'"'"-'(a;, x'; t) is expressed as 

/ + OC 
dk F{x,x']k)e-"'^\ (29) 

with Tf defined by 

ht 

Tt = T^. (30) 

Here, x'; k) is defined by 



2M { 1 - cos[A:(x + a;' - 2L)] + 2f sin[fc(a; + x' - 2L)] ■ 
F{x,x'-k) = -^<^cos[fc(x-.T')] + ^ ^ ^ ^ \ (31) 



1 4- ^^ 



with 



C{k)^ 



A 

'k 



(32) 



It may be noted that the integral region of the wave number 
k in the right-hand side of Eq. ( |29] l is (— oo, +oo), different 
from in the right-hand sides of Eqs. (|26] |. 

Noting that the function ( 1311 1 is an even function of k, we 
expand it as 



+ CJO 



F{x,x';k) = ^B^{x,x')k^'' 



(33) 



u=0 



with respect to fc^, where B^{x, x') is a function of x and x' 
and is independent of k. Then, for t > the fc-integral in Eq. 
( |29] l can be carried out, and we obtain 



G(«^°)(x,x';0 = ^ 



noting 0! = 1. Here, we used the integral formula 



2u+l ■ 



dk k 



2v —tTth 



{2v)\\-i j — ¥ 



(34) 



(35) 



for V = 0, 1, • • • , which is derived from the rt-derivative of 
the equation dk e'"*^' = (1 - i)^/Tr/{2Tt) for t > 0. 
Eq. ( |34] | is an asymptotic expansion of the Green function 
G'(°'°)(x,x';i) with respect to l/t'^+i/^, = 0, 1, - • •. 



Under the initial condition (|27] | we describe escape behav- 
iors of a particle from the finite region As an example of 
quantities to characterize such escape behaviors, we consider 
the quantity P{t) defined by 



^ J,'dxp(^\x,t) 
/>xp(o)(x,0)' 



(36) 



This is the ratio for a particle to survive in the finite region 

£(0) 

at time t in comparison with that at the initial time t ~ 
1I47I1 . and we call this probability P(t) the survival probability 
in this paper ll48ll . As another quantity to characterize particle 
escape behaviors we also consider the local velocity for the 
particle to escape from the finite region C^^\ which is defined 
by 



V{t) 



= - lim u(°)(x,t) 



(37) 



with the probabiHty density p{0,t) = |^'(0,i)p at the junc- 
tion O. Here, we used Eq. (|5]l to derive the second equation 
in the right-hand side of Eq. ( |37] ) from its first equation. We 
call this velocity V{t) the "escape velocity" in this paper The 
survival probability ( |36] | and the escape velocity ( |37] | can be 
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calculated by the wave function ( |28] | only in the finite region where C{k) is given by 
£(0) via the Green function ( l34b . without information on the 
particle in the semi-infinite region 

In the following subsection, we consider properties of the 
survival probability P{t) and the escape velocity V{t) mainly 
in the fixed boundary case /i — > +00, while we give some 
analytical arguments of P{t) and V{t) in the open boundary 
case /i = in AppendixiDl 



B. Asymptotic properties of particle escapes 

In the fixed boundary case yu +00, the function ( |3TI ) is 
represented as 

x^fc)^!^ (38) 



C(fc) 



TV + 1 - (TV - 1) cos(2fcL) + A 



sin(2fcL) 



k 



{N - 1) siii(2fcL) - A 



1 - cos(2fcL) 



(39) 



The expansion of the function ( |38] | with respect to as in Eq. 
(|33] | is different between the cases of AL 7^ —1 and AL = — 1, 
so we consider these two cases separately below. 



J 



First, we consider the case of XL 7^ —1. For the function 
F{x, x'\ k) with respect to k"^^ in Eq. ( |33] | is given by 

Bo{x,x') = 0, 

N{x-L){x' -L) 



Bi{x,x') 
B2{x,x') 



7r(l + AL)2 ' 
N{x^ L){x' - L) 
67r(l + Xiy 



1. Fixed boundary case with XL 7^—1 



the quantity B^{x,x') as the coefficients of the function 

(40) 
(41) 



(x - Lf + {x' - Lf + 2L 



^?,{N^ - I) - \L{\L + A) 
(1 + AL)2 



(42) 



• • • , concretely. By inserting these coefficients Bi,{x, x'), = 0, 1, • • • into the formula (|34l l and by using Eq. ( |28T l we obtain 
the wave function \E'('^' {x, t) in for i > as 



^^°'>{x,t) 



N{x~ L){l + i)ei 1 



3/2 



2V2^(1 + AL)2 
A^(a;-L)(l-i) 1 
' 8V2^(1 + AL)2 /2 



(43) 



where is defined by By Eqs. ([T]i, ( l30b and ( |43] | the survival probability ( l36l l in 

2, the long time limit is expressed asymptotically as 

e^= dx {x- L)''¥°\x,0). (44) 
^0 

Eq. (03]) is the wave function for the fixed boundary case 

fj, +00 with XL =/= —1 as an expansion of l/t"^^^'^, , , 2N'^Ei / mi\ 

i^=l,2,---. ^^^^ ^ 37r(l + AL)4 [l^J ^^^^ 
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with the constant Si = 8f / /J^ dx p^'^^x, 0) determined by 
the initial wave function ^'("^(x, 0) only. Eq. ( |45] l means that 
the survival probability P{t) decays in power ^ i"'^ asymp- 
totically in time t. Moreover, by Eq. ( |45] | we obtain 



Pit) 



lim - 



(46) 
(47) 



^ of has already been 
1 in Ref. lElll. On the other 



as far as the initial wave function ^'^"^ {x, 0) is independent of 
A and A^. Eq. ( |46] ) shows that the probability for a particle 
to stays in the finite region becomes lower for a larger 
amplitude (|7]i of potential barrier at the junction O in the long 
time limit, and this (rather counter-intuitive) result, as well as 
the asymptotic power decay ^ t 
shown for the one-lead case N ~ 
hand, Eq. ( l47l ) means that the probability for a particle to stays 
in the finite region £(0) for multi semi-infinite leads becomes 
A^^ times higher than that for a single lead. 

The first non-zero contribution of the survival probability 
P{t) after a long time comes from the first term of the right- 
hand side of Eq. ( |43] ). In contrast, the escape velocity V{t) 
is an example of quantities in which the first non-zero contri- 
bution after a long time comes from the second term of the 



right-hand side of Eq. ( l43T l. Actually, by inserting the right- 
hand side of Eq. (03]) up to its second term into Eq. (O, and 
using Eq. ( l30b we obtain 



v^°\x,t) 



t— J- + CX) 



L ~ X 



(48) 



asymptotically in time. Therefore, after a long time the escape 
velocity dSTI l is represented as 



V{t) 



L 



(49) 



meaning that asymptotically in time, the escape velocity V{t) 
decays in power ^ t^^ and is independent of the number of 
semi-infinite leads, the constant A, and the initial wave func- 
tion \I'("'(a;, 0), etc. The velocity (|49] l can be regarded as a 
constant velocity by which a classical mechanical particle in 
an ideal wire remains inside a finite region with the length L 
within the time interval t without an escape. 



2. Fixed boundary case with XL = — 1 

Next, we consider the case of XL = — 1, in which the quan- 
tities Bt,{x, x'), ly = 0,1, ■ ■ ■ are represented as 



Bo{x,x') 
Bi{x,x') 



{x-L){x' -L) 
{x-L){x' - L) 

mTl? 



{x - Lf + {x' - Lf 



2L2 37v2 _ 1 



iV2 



(50) 
(51) 



• • • , concretely, by Eqs. (|33T l, ( [38] l and (|39T l. By inserting these coefficients B^{x, x'), v ~ - ■ ■ of the function F{x, x'; t) 
with respect to fc^ into Eq. (|34] | the wave function ( |28] | in the region for t > is represented as 



{x - L){1 - i)<di 1 {x~L){l + i) 



2TrNL^ 



1/2 



12V2^NL^ 



{x - LfQi + 63 



2L^ 7V2 _ 1 



iV2 



3/2 



(52) 



as an expansion with respect to l/V^^^/"^, = 0, 1, • • • . 

From Eqs. ([T]i, ( l30l ). (|36] | and ( |52] | we derive the survival 
probability 



P{t) 



SnhN^Lt 



(53) 



asymptotically in time. It is important to note that the survival 
probability (|53] | for XL = —1 decays in power ^ t~-^ asymp- 
totically in time, different from the asymptotic decay power 
~ for XL 7^ —1 as shown in Eq. ( |45] ). By the asymptotic 



(54) 



form of the survival probability (1531 ) we obtain 



lim — = — - 



for the A^-independent initial wave function ^'("^(a::, 0). Eq. 
( [54b means that by connecting more number A^ of the leads to 
the finite region C^'^K the probability of a particle to stay in 
the region C^^"* for XL ~ —1 becomes lower after a long time, 
opposite to Eq. (|47] | for XL 7^ —1. 

From Eqs. (O, (|30] | and ( |52] | we derive the local velocity 
in£(o) as 



L ~ X 

it 



(55) 



asymptotically in time. Therefore, we obtain the escape ve- 
locity dJTl i as 



V{t) 



L 
3I 



(56) 
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after a long time. Eq. ( |56] ) shows that the escape velocity V{t) 
for XL = — 1 decays asymptotically in the same power ^ t^^ 
as in Eq. ( |49] l for XL 7^—1, although its prefactor i/3 is one 
third of that in the case of XL 7^—1. 

C. Finite time properties of particle escapes 

In this subsection, we consider a finite time behavior of par- 
ticle escapes in the fixed boundary case /i +oo by calcu- 
lating the wave function 'i>^'^\x,t) numerically. In order to 
calculate the wave function 5''^"^ (x, t) concretely for the fixed 
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FIG. 3: Survival probability P(t) as a function of t for the fixed 
boundary case with XL 7^—1 for N = 1 (the solid line), 2 (the 
broken bne), 4 (the dotted line), 8 (the dash-dotted line) and 16 
(the dash-double-dotted line). The main figure is log-log plots of 
P{t) with the corresponding asymptotic power decays shown in thin 
straight lines, and the inset is their linear-log plots for a short time. 



with C{k) given by Eq. ( [39] ). By carrying out the integral ( fSST l 
with respect to k numerically, we calculate the serval proba- 
bility P{t) and the escape velocity V{t) for a finite time t. 
In this subsection we chose the parameter values as L = 1, 
tr = 1, TO = 1 and h — 1. 



boundary case, we specify the initial wave function in the fi- 
nite region 

£(0) 

as 

*^°^(a;,0) = y|sin(^x) , (57) 

which is the cr-th eigenstate of a particle confined in the finite 
region £(0) without semi-infinite leads (a = 1, 2, • • • ) 114911 . 
Under the initial condition (ET]) and (US, by Eqs. (|28l), (|29] l 
and (l38l l the wave function vI'C') (x, t) in the finite region C^^^ 
at time t is represented as 



(58) 

I 
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FIG. 4: Absolute value of the escape velocity as a function 

of t for the fixed boundary case with XL 7^—1 for A'^ = 1 (the solid 
line), 2 (the broken line), 4 (the dotted line) and 8 (the dash-dotted 
line). The main figure is log-log plots of with thin straight 

lines to show the asymptotic power decays, and the inset is linear- 
linear plots of the escape velocity V(t) itself for a short time. 



L Fixed boundary case with XL 7^—1 

Figure |3] is the survival probability P{t) as a function of 
time t for the fixed boundary case fj, +oo with XL 7^ — 1 
under the initial wave function ( fSTl i for different numbers 
iV = 1 (the soHd line), 2 (the broken Une), 4 (the dotted 
line), 8 (the dash-dotted line) and 16 (the dash-double-dotted 
line) of semi-infinite leads. Here, we used the parameter value 
A = 1, and the thin straight lines in this figure show the 
asymptotic power decay ( |45] | of the survival probability for 



^—irt k 



dk 



sin(fcL) sin[fc(a; — L)] 
C{k) [{kLY - (cr7r)2] 
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each value of N. 

The numerical results in the main figure of Fig. [3] as log- 
log plots of P{t) show that the survival probability P{t) ap- 
proaches to the power decay ( |45] | asymptotically in time from 
an earlier time for a smaller number N of leads. We can also 
see in the inset of Fig. |3]as hnear-log plots of P{t) that for a 
short time the survival probability P{t) decays exponentially 
^ cxp{—at) in time with a positive constant a, as shown es- 
pecially in the cases of = 4, 8 and 16. The time period for 
such an exponential decay of P{t) is longer for more number 
iV of leads, and its decay rate a is smaller for more number N 
of leads. One may notice that the survival probability does not 
decrease monotonously in time and shows a time-oscillation 
between its exponential decay region and the power decay re- 
gion. As a tendency, the survival probability P{t) is higher 
for more number N of leads, although there are temporal ex- 
ceptions for it by its time-oscillatory behavior. 

Figure|4]is plots of the absolute value \V{t)\ of the escape 
velocity as the main figure, as well as the escape velocity V{t) 
itself as the inset, as a function of time t for the fixed bound- 
ary case under the initial wave function (|57] | for the different 
lead number = 1 (the solid line), 2 (the broken line), 4 
(the dotted line) and 8 (the dash-dotted line). Here, we used 
the parameter value A = 1, and the thin straight line is the 
asymptotic form ( |49] l of the escape velocity. 

It is shown in the inset of Fig. |4]that in a short time the 
escape velocity V{t) oscillates in time, rather than a simple 
decay, and can take even a negative value sometimes. Such a 
time-oscillatory behavior of V{t) continues for a longer time 
for more number N of leads, although its time oscillating pe- 
riod seems to be almost independent of N. For a short time, 
as a tendency the magnitude of the escape velocity V{t) be- 
comes larger for more number N of leads. For some values of 
N, such as for = 1 in the inset of Fig. |4] a very rapid (but 
not abrupt) change of the escape velocity V{t) from a posi- 
tive value to the first negative value as a function of t occurs, 
when the value of |v]>(") (+0, t)\ is non-zero but very small at 
the time t satisfying the condition 94'^"-'(x,i)/9a;|x=+o = 0. 
(Note that time-oscillations of V{t) around the time i « 10 in 
the main figure of Fig. |4]look like to be cut off in a middle 
occasionally, because there are not enough numbers of cal- 
culation points, but they actually reach the value of zero by 
crossing the line V = 0.) After such a time-oscillation, the 
escape velocity V{t) converges rapidly to its asymptotic form 
(|49] |. which is independent of N. 

2. Fixed boundary case with XL = — 1 

Now, we consider finite time properties of particle escapes 
for XL = —1 (so A = —1 because of L = 1) in the fixed 
boundary case /^t +oo. Figure|5]is the survival probability 
P{t) as a function of time t for = 1 (the solid line), 2 (the 
broken line), 4 (the dotted line), 8 (the dash-dotted line) and 
16 (the dash-double-dotted line). The thin straight lines in this 
figure show the asymptotic form i53[ of P{t) for each value 
ofN. 
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FIG. 5: Survival probability P{t) as a function of t for the fixed 
boundary case with XL = — 1 for N — 1 (the solid line), 2 (the 
broken line), 4 (the dotted line), 8 (the dash-dotted line) and 16 
(the dash-double-dotted line). The main figure is log-log plots of 
P{t) with the corresponding asymptotic power decays shown in thin 
straight lines, and the inset is their linear-log plots for a short time. 




FIG. 6: Absolute value |'t^(t)| of the escape velocity as a function 
of time t for the fixed boundary case with XL = — 1 for A'^ = 1 (the 
solid line), 2 (the broken line), 4 (the dotted line) and 8 (the dash- 
dotted line). The main figure is log-log plots of with a straight 
line for the asymptotic power decay, and the inset is liner-linear plots 
of the escape velocity V(t) itself for a short time. 



As shown in the inset of Fig. |5]as linear-log plots of P{t), 
especially for large A^ = 8 and 16, in a short time region the 
survival probability P{t) decays exponentially in time. Such 
an exponential decay continues for a longer time for more 
number A^ of attached leads, and its decay rate is smaller for 
more number A^ of leads. After the exponential decay, the 
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survival probability P{t) shows a time-oscillatory behavior, 
then converges to its asymptotic power decay form ( l53b , as 
shown in the main figure of Fig. |5]as log-log plots of P{t). 
Contrast to the case of XL ^ —1, the survival probability 
P{t) becomes lower for more number of N in this asymptotic 
power decay region, while it is opposite in the exponential 
time-decay region. 

Figure|6]is the absolute value \ V{t) \ of the escape velocity 
as the main figure, as well as the escape velocity V{t) itself 
as the inset, as a function of time t for the fixed boundary 
condition with XL = — 1 in the case of = 1 (the solid line), 
2 (the broken Une), 4 (the dotted line) and 8 (the dash-dotted 
line). The thin straight line in this figure is the asymptotic 
form ( |56] | of the escape velocity V{t). The escape velocity 
V{t) itself is also shown for a short time in the inset of this 
figure. 

As shown in the inset of Fig. |6l the escape velocity V{t) 
oscillates in time, first as a time-oscillation with positive val- 
ues then as that with positive and negative values. Such an 
oscillatory time region continues for a longer time for more 
number N of leads, although its time oscillating period seems 
to be almost independent of A^, and as a tendency the absolute 
value \ V{t) \ of the escape velocity is larger for more number 
N of leads in this time region. (In the main figure of Fig. |6] 
the time-oscillations of V{t) look to be cut off in a middle, 
because of little number of calculation points. In the actual 
graphs the value of \V{t)\ go to zero in a time-oscillatory re- 
gion with positive and negative values.) After such a time- 
oscillation the escape velocity V{t) converges rapidly to its 
asymptotic power decay form ( |56] |. which is independent of 
N. 



IV. CONCLUSION AND REMARKS 

In this paper, we have discussed particle escapes in an open 
quantum network system by a scattering theoretical approach. 
As a concrete example with a current splitter as a feature of 
network systems, we considered particle escapes from an end 
of a finite one-dimensional wire to N number of semi-infinite 
one-dimensional leads. Properties of particle escapes in such 
a quantum network was discussed by using the two kinds 
of quantities; the one is the probability P{t) for the particle 
to remain in the finite wire at time t, the so-called survival 
probability, and the other is the velocity V{t) for the parti- 
cle to leave from the finite region, the so-called escape ve- 
locity. Here, the escape velocity V{t) is introduced from the 
probability current, based on the equation of continuity for the 
particle position probability density. With the fixed boundary 
condition at an end of the finite lead, for the potential bar- 
rier amplitude A 7^ — 'h?/{2mL) the survival probability P{t) 
depends on the number N of attached semi-infinite leads as 
\\mt-^j^ryo P{t) / P{t)\N=i = and decays in power ^ 
asymptotically in time. In contrast, for the potential barrier 
amplitude A = —h^/{2mL) the survival probability satisfies 
the relation \im.t^+oo P{t) / P{t)\N=i = and it decays 

in power ^ \/t after a long time. On the other hand, the es- 



cape velocity V{t) decays like CL/t asymptotically in time 
with the constant C which is independent of the number N 
of leads and the initial wave function ^''^^(a;, 0), etc. It was 
also shown that for a finite time the survival probability P{t) 
decays exponentially in time for a longer time with a smaller 
decay rate for more number N of attached leads, and shows a 
time-oscillatory behavior between the exponential decay time 
region and the power decay time region. The escape velocity 
V{t) show a time-oscillatory behavior for a finite time, and as 
a tendency its value is higher for more number N of attached 
leads with a larger amplitude of time-oscillations. 

We described the dynamics of an escaping particle by a 
quantum scattering theoretical approach. It may be noted 
that although the time-dependent wave function of an escap- 
ing particle is expanded by the scattering states in the open 
network system it is nonzero only in a finite region for a fi- 
nite time and is normalizable, different from stationary quan- 
tum scattering states caused by incident plain waves from in- 
finitely far spatial regions. In order to construct concretely the 
scattering states in the open network system, it is essential to 
specify the boundary conditions at the junction of a finite wire 
and multiple leads and at another end of the finite wire. We 
specified these boundary conditions based on the probability 
current given from the equation of continuity for the particle 
position probability density. We imposed the conservation of 
this current and the continuities of the particle wave function 
at the junction, so that the scattering matrix at the junction is 
unitary and the scattering states satisfies the orthogonal rela- 
tion automatically. It is important to note that the condition 
of no net current at an end of the finite wire leads to a group 
of boundary conditions specified by a parameter ^. In this 
paper we mainly considered the case of /i — > +00, i.e. of 
the fixed boundary condition. However, by using the case 
of /z = 0, i.e. of the open boundary condition, we can also 
discuss particle escapes from a finite one-dimensional wire 
whose both ends are connected with multiple one-dimensional 
semi-infinite leads. In this case, for the initial state which 
is anti-symmetric with the reflection at the center of the fi- 
nite region, we get the same results of the survival probabil- 
ity P{t) and the escape velocity V{t) as cases of the fixed 
boundary condition. In contrast, for the initial state which 
is symmetric with the reflection at the center of the finite re- 
gion, by using the open boundary condition /i = we obtain 
limt^+ooP(i)/[P(t)U=i] = and 1/(0 3L/t 
for A ^ 0, and Vmit^+oo P{t) /[P{t)\N=i] = l/N^ and 
V{t) Lit for A = 0, in which different behaviors of 

P{t) and V{t) occur at a different value of A from that in the 
anti-symmetric initial state. As a remark it may be interesting 
if we could clarify the physical meanings of particle escapes 
in the case of other values of the parameter /i, i.e. a nonzero 
and finite values of /i. We also note that the current used to 
specify these boundary conditions is a current of the probabil- 
ity density of the particle position. In this sense, the escape 
velocity V{t) based on this probability current at the junction 
is not the particle velocity itself. The uncertain principle of 
quantum mechanics forbids to specify the particle velocity at 
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FIG. 7: Escape velocity V{t) as a function of t for the fixed bound- 
ary case with A = (the solid line), 4 (the broken line), 8 (the dotted 
line) and 12 (the dash-dotted line). 
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Appendix A: Scattering matrix 

In this appendix we give a derivation of Eq. ( fTTT i with the 
scattering matrix (fT2b . 

By inserting Eq. ( fTol i into the condition ^ we obtain 



^'(0,<). (A.l) 



On the other hand, Eqs. (|6]l and (fTOl) lead to 



a specific position of particle. 

As finite time properties of quantum escapes, in this paper 
we discussed mainly dependences of the survival probability 
P{t) and the escape velocity V{i) on the number N of at- 
tached leads, but their dependences on other parameters also 
express some important escape properties. For example, the 
survival probability P{t) decays exponentially for a finite time 
with a smaller decay rate for a larger value of the parameter 
A proportional to the potential barrier amplitude, as known al- 
ready for the case of iV = 1 ll22i I25I1 . As another example, 
we show Fig. |2]for the escape velocity V{t) as a function of 
time t for the fixed boundary case /i +00 with A = 0, 4, 8 
and 12. Here, we chose the other parameter values as L = 1, 
= 4, (T = 1, 771 = 1 and h = 1. This figure shows that the 
escape velocity V{t) oscillates in time around a constant value 
Vq for a short time, and the value Vq is almost independent of 
A but the time period appearing such time-oscillations around 
the value Vq becomes longer for a larger value of the param- 
eter A. Detailed dependences of escape behaviors in quantum 
networks via multiple leads on other parameters, such as A, fi 
and (7, etc., will be discussed elsewhere. 

In this paper we discussed briefly exponential decays of the 
survival probability for a finite time, appearing especially for 
a large number of attached leads as well as for a large value 
of the potential barrier amplitude A. These exponential decays 
would be related to the poles of the Green function or the scat- 
tering matrix i3l ll9ll2ll 122113511 . A detailed analysis of such 
exponential decays of the survival probability as the scattering 
resonance in particle escapes via multiple leads remains as an 
important future problem. 



N 



) _^(n) 
k) °(i,fc) 



n=0 



(A.2) 



By Eqs. 



A. lb and (IA.2ll we obtain 



2a 



(n) 



2 



N 



*(0,i) 



(n) 



(A.3) 



From Eqs. dA.U and ( IA.3I ) we derive 

= *(0,i)(ll...l)^ 



(A.4) 



with the scattering matrix ( [12] ). Therefore, we obtain Eq. ( fTTI ). 



Appendix B: Ortliogonality of scattering states 



In this appendix we give a proof of the orthogonahty rela- 
tion ( I2TI ) for the scattering state ^''(^\) (x)- 
We note the mathematical identity 



r + oo 



+00 

dxe'^-" = lim / dxe**^^-'^ 



7r(5(fc) + ix{k) 



where S{k) is the delta function 
S{k) = lim — - 



»+0 TT fc2 _|_ £2 



(B.l) 



(B.2) 
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and x(fc) is defined by 



where we used Eq. (IrTT i and the identity (1 - e^^^)8{k) = 0. 



X{k) = lim 



k 



Using the function ( IB. 31 ) we obtain 



(1 - 6*'=^) ix{k), 



(B.3) 



By the scattering state ( fTOl i and the inner product (|20] | for 
quantum states, as well as Eqs. (fTTT i. (fT4l i. ( IB. II ) and ( IB.4l i. we 
(B.4) obtain 



N 



n=l •^'^ 

N 

<5(fc-fc')+-EH;i)*^K.')+^ 



"(i^fe) '^OAfe') ^°{j,k) °(o',k') 



Tl=l 



e->+0 
TV 

n=l 

- /) x{k - k') + (5^ - Su) x{k + k') 

r 



(n) * (n) 



(0) 



(0) ^-i(k-k')L 



S{k + k') 



(B.5) 



in which is defined by 



Sk 



(B.6) 



Now, we note 

AT 



n=l 



(n) * (") 



5{k-k') 



= 2 



„„(0) * (0) _,(0) *,(0) 1 



TV 



n=l 



fc') 



(B.7) 



by Eqs. dTT), dlSl, (O, (HSll and (|22li. Moreover, using the 
function ( IB. 3b we note 



& x{k)X{k) =V dk 



X{k) 



for any function X(fc) of fc, where we introduced the operator 
V as that to take the principal integral. In this sense, for the 
scattering matrix (fT2l) we obtain 



(slSk- - l) x{k - k') + (5^ - Su) x{k + k') 



= O 



SlSk' - 1 



k - k' 



k + k' 



(B.9) 



with the matrix O whose all elements are zero. Similarly, for 
the quantity Sk defined by Eq. iB.6i we obtain 



{slsk'-l)x{k-k') + isl-sk')x{k + k') 
= V 
= 0. 



(B.IO) 



(B.8) Inserting Eqs. dBJ]) . dBlgj l and dfilOt into Eq. (Ib3] i we obtain 



TV 



(n) *, (n) , (n) * (n) 



5{k + k'). 



(B.ll) 
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Noting that the second term of the right-hand side of Eq. By Eqs. ( fTOb and ( fT4b we note 

( IB.l lb is zero for kk' > 0, we obtain the orthogonal relation 
( l2n i of the scattering states. 



Appendix C: Green function G^"'"' {x, x'; t) 



In this appendix we show Eq. ( |29] ). 





2 







1+2^ 

ik{x — L) _|_ —ik(x — L) 



= 2 

4 

TT 



,(0) 



cos[fc(a; — x')] + Re 



I - it 

^ — ik{x — L) _j_ k_^ik{x —L) 



1 + ii 



1 -i 



k_^ik{x+x — 2L) 



cos[A;(x — a;')] + 



cos[fe(3;+x'-2L)]+2-ii sin[fc(2;+a;'-2L)] 



iV + 1 + if + (TV - 1 + j|) i±^e2i''i 



(C.l) 



where we used Eqs. ( [TtT i. ( fT9] l and ( |22] |. Here, Re[X] is the 
real part of X for any complex number X. It is important to 
note that the quantity (IC.lb is an even function of k, namely 

'5'(^.)(-)'5'(^.)(-T - $g),)(.)$g;,)(xr. (C.2) 

By inserting Eq. (IC.lb into Eq. (l26T l for n = n' ^ {) and by 
noting Eq. (|C2] | we obtain Eq. (|29ll with Eqs. dlOll and (TO . 



In this case, the expansion of the function ( ID. lb with respect 
to as in Eq. (|34] | is different between the cases of A = 
and A 7^ 0, so we discuss these two cases separately below. 
1. Open boundary case with A / 



Appendix D: Escape beliaviors in the open boundary case 



In this appendix we discuss asymptotic escape behavior of 
a particle in the open boundary case = 0. In this case, the 
function ( 13Tb is represented as 



F{x,x';k) 



4iV cos[fc(a: - L)] cos[fc(a;' - L)] 



C{k) 



(D.l) 



where C(k) is given by 



C{k) = 



N + 1 + {N -I) cos{2kL) - A 



sin(2A:L) 



{N - l)sin(2fcL) + A 



1 + cos{2kL) 
k 



(D.2) 



In the open boundary case /i = with nonzero potential 
barrier with A 7^ at the junction O, the expansion coef- 
ficients B^{x,x'), V ~ 0, 1, • • • of the function (ID. lb with 
respect to fc^, like in Eq. (l33b . are concretely represented as 
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Bo{x,x') = 0, 
Bi{x,x) = 



B2{x,x') = - 



N 
27rA2 



/ , r^2 - 2AL- (AL)2 



A2 



• • • , concretely. By the Green function (|34] | with these coefficients B^{x, x'), = 0, 1, • • • , the wave function 
expressed as 



iV(l + z)eo 1 3jV(l-i) 



7V2 - 2Ai: - ( AL)2 ■ 
(x - Lfe^ + 92 + 2 ^ ^ 9o 



1 



5/2 



r 



(D.3) 
(D.4) 

(D.5) 

for t > is 

• • • (D.6) 



with &^ is defined by Eq. ( |44] |. 

By Eq. ( ID.6I 1 the survival probability P{t) is represented 
asymptotically in time as 



P{t) 



(D.7) 



7rA4 

with the constant Sq = Og/ /^^ c?x p'-°^(a;, 0) determined by 
the initial wave function \[''^°)(a;, 0) only. Eq. (ID.7b leads to 
the relations 



lim , , , 
lim — -— ^ 

t->+oo P(t)|A=l 



1 

A4 



(D.8) 
(D.9) 



for the N- and A- independent initial wave function 
(x, 0). Therefore, the A^-dependence ( ID. 81 ) of the asymp- 
totic survival probability for the open boundary case with 
A ^ is the same as that shown in Eq. WT\ for \L ^ —1 
in the fixed boundary case. In contrast, the asymptotic sur- 
vival probability is inversely proportional to A"*^ in Eq. ( ID. 71 ). 
while in the fixed boundary case with XL 7^ —1 it is inversely 
proportional to (1 + AL)* as shown in Eq. ( |45] |. 



From Eqs. (O and ( ID.6b we also derive an asymptotic ex- 
pression of the local velocity f (x, t) in the finite region 

£(0) 

as 



3(L-a;) 



t 



(D.IO) 



Therefore, the escape velocity (l37b is given by 

t-S-+oo 3L 



t 



(D.ll) 



asymptotically in time, meaning that the escape velocity de- 
cays in power ~ as in the fixed boundary case, but its 
value is three times as high as in Eq. ( |49] l. 



2. Open boundary case with A = 

In the case of zero potential barrier with A = 0, from 
Eqs. (l33T l and (ID. lb we derive the quantities By{x,x'), v = 
0, 1, • • • as 



Biix.x') - 



2ttN 



{x - Lf + {x' - Lf - 2L 



.N^ - 1 

iV2 



(D.12) 
(D.13) 



• • • , concretely, in the open boundary case /i = 0. By using these coefficients Bi,{x, x'), = 0, 1, • • • and Eqs. ( l28T l and ( |34] |. 
the wave function ^'^"^ (cc, <) for i > in is represented as 



(l-i)Bo 1 1 + i 



2iiN tI'"^ 4V2^iV 



(x - LfQo + 62 - 2L 



2^'-l ■ 



3/2 



(D.14) 



r 



with defined by Eq. ( |44] |. 



By Eq. ( ID. 141 ) the survival probability P{t) is represented 
asymptotically in time as 

t->+oo 2^QmL 



Pit) 



ttN^M 



(DAS) 
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Eq. ( ID.I5I 1 shows that the survival probability P{t) in the case 
of /^t = and A = decays asymptotically in power ~ tr^, 
as in Eq. ( |53] | for the fixed boundary case with XL = — 1. 
Moreover, from Eq. ( ID.15b we derive the relation with the 
same as Eq. ( l54l i for the A and A^-independent initial wave 
function \['(°)(x, 0). 



From Eqs. (|3]l, (|30] | and ( ID.14l i we derive the local velocity 
v''^^ (x, t) in the finite region which is asymptotically the 
same as Eq. (|48] | in time. Therefore, we obtain the same 
asymptotic escape probability V{t) as Eq. ( |49] l for the fixed 
boundary case with AL 7^—1. 
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